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We investigate transport properties of a superconducting junction of many (A'^ > 2) one- 
dimensional quantum wires. We include the effect of electron-electron interaction within the one- 
dimensional quantum wire using a weak interaction renormalization group procedure. Due to the 
proximity effect, transport across the junction occurs via direct tunneling as well as via the crossed 
Andreev channel. We find that the fixed point structure of this system is far more rich than the 
fixed point structure of a normal metal— superconductor junction {N = 1), where we only have 
two fixed points - the fully insulating fixed point or the Andreev fixed point. Even a two wire 
(A = 2) system with a superconducting junction i.e. a normal metal— superconductor— normal 
metal structure, has non-trivial fixed points with intermediate transmissions and reflections. We 
also include electron-electron interaction induced back-scattering in the quantum wires in our study 
and hence obtain non-Luttinger liquid behaviour. It is interesting to note that (a,) effects due to 
inclusion of electron-electron interaction induced back-scattering in the wire, and (h) competition 
between the charge transport via the electron and hole channels across the junction, give rise to a 
non-monotonic behavior of conductance as a function of temperature. We also flnd that transport 
across the junction depends on two independent interaction parameters. The first one is due to the 
usual correlations coming from Friedel oscillations for spin-full electrons giving rise to the well-known 
interaction parameter {a = ((72 — 2g\) /2TihvF)- The second one arises due to the scattering induced 
by the proximity of the superconductor and is given by (a' = ((?2 + gi) /2nhvF). The non-monotonic 
conductance and the identification of this new interaction parameter are two of our main results. In 
both the expressions gi = V{2kF), g2 = 1^(0), where V{k) is the inter electron interaction potential. 

PACS numbers: 73.23.-b,74.45.+c,71.10.Pm 



I. INTRODUCTION 

Effects due to the proximity of a superconductor has 
motivated a lot of worbii^i^ in the last several decades. 
A direct manifestation of proximity effect is the phe- 
nomenon of Andreev reflection (AR) in which an electron 
like quasi-particle incident on normal— superconductor 
(NS) junction is reflected back as a hole along with the 
transfer of two electrons into the superconductor as a 
Cooper pair. An even more intriguing example where 
the proximity effect manifests itself is the phenomenon 
of crossed Andreev reflection (CAR) which can only take 
place in a normal metal— superconductor— normal metal 
(NSN) junction provided the distance between the two 
normal metals is less than or equal to the phase co- 
herence length of the superconductor. This is a non- 
local process where an incident electron from one of the 
normal leads pairs up with an electron from the other 
lead to form a Cooper pair— i^i^iii^ and joins the super- 
conductor. Its relevance in the manipulation of spin 
currents^. (SC) and questions regarding production of 
entangled electron pairs in nano devices for quantum 
computation has attracted a lot of attention in recent 
times^iiii^ii^i^ii^i^ii^i^. Further extensions such as 
inclusion of effects due to electron-electron interactions 



on AR processes in case of NS junctions in the context of 
one-dimensional (1-D) wires have also been considered 
recentlj*i2^^i22^'24. 

In this article, we shall first develop a general formu- 
lation for studying the transport properties of a mul- 
tiple quantum wire (QW) junction in the spirit of the 
"Landauer-Buttiker" approacb^^, where the junction it- 
self is superconducting. We shall use this formulation to 
study the influence of the proximity effect on the trans- 
port properties of a superconducting junction speciflcally 
for the case of two and three 1-D interacting quantum 
wires and show how the simple case of junction of a single 
1-D QW with a superconductor (NS junction) is different 
from the multiple wire junction counterpart. Because of 
the existence of the AR process, both electron and hole 
channels take part in transport. The power law depen- 
dence of the Andreev conductance for the NS junction 
case was first obtained using weak interaction renormal- 
ization grou p (WIRG) approach by Takane and Koyama 
in Ref. Il9l . This was in agreement with earlier re- 
sults from bosonization™, which, however, could only 
handle perturbative analyses around the strong back- 
scattering (SBS) and weak back-scattering (WBS) limits. 
The WIRG approach, on the other hand, can study the 
full cross-over from WBS limit to SBS limit. Hence the 
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FIG. 1: Multiple wires connected to a superconducting junc- 
tion. The dashed lines represents the fact that the model 
can be trivially extended to more than two wires, 'a' is the 
effective length of the superconductor. 



WIRG approach is very well-suited for studying problems 
where the aim is to look for non trivial fixed points with 
intermediate transmissions and reflections. This would 
be difficult using a bosonization approach. 

For the NS junction case, it was shown that the power 
law exponent for the temperature dependence of conduc- 
tance was twice as large as the exponent for a single bar- 
rier in a QW. This happens because of the introduction 
of the extra hole channel due to AR. The WIRG approach 
takes into account both (a) electron-electron interaction 
induced forward scattering processes which gives stan- 
dard Luttinger liquid behavior, and (b) electron-electron 
interaction induced back-scattering processes which give 
rise to non Luttinger liquid behavior (deviation from pure 
power law behavior) as obtained by Matveev et al. in 
Refs. I26I and (2^. The most interesting point to be noticed 
here is the fact that both bosonization and WIRG give the 
same Luttinger liquid power law dependence for the con- 
ductance of NS junction although the WIRG approach 
actually takes into account the extra process of electron- 
electron back-scattering which usually leads to non Lut- 
tinger liquid behavior. This happens because there is 
a remarkable cancellation in perturbation theory which 
nullifies any deviation from pure power law behavior of 
conductance. 

In this article, we show that this kind of cancella- 
tion does not happen for the NSN junction or for that 
matter, for any junction comprising of more than two 
1-D quantum wires. Hence, deviations from pure power 
law do exist. We show that due to the inter-play of the 
proximity and the interaction effects, one gets a novel 
non-monotonic behavior of conductance for the case of 
NSN junction as a function of the temperature. Note that 



this is something which cannot be obtained in a bosoniza- 
tion analysis which neglects the back-scattering part of the 
electron- electron interaction. We extend our results to 
junctions with a ferromagnetic wire on one side, i.e., fer- 
romagnetsupcrconductor-normal (FSN) junctions and 
to junctions with fcrromagnets on both sides, i.e., ferro- 
magnet-superconductor- ferromagnet (FSF) junctions. 
Here we assume that the influence of the junction be- 
tween the superconductor and the ferromagnetic wire has 
a very small effect on the spectrum of the superconductor 
itself. Of course this is true only if the superconductor 
is large enough. We also study transport through a su- 
perconductor at the junction of three wires (and finally 
extend it to N wires). This generalizes the earlier work 
on junction o^^'^^ to also include proximity effects. 

In Sec. |TT1 we review the applicability and strength of 
the WIRG approach when applied to quantum impurity 
problems, such normal junction of multiple quan- 

tum wires, or a junction with a spin impurity. Then we 
discuss how one can apply the WIRG technique to the 
problem studied in this paper. In Sec. Illli we describe 
the set-up for our system, i.e. a superconductor at the 
junction of N wires, in terms of a scattering matrix and 
briefly discuss the symmetries of the proposed model. 
We then perturbatively (in electron-electron interaction 
strength) calculate the leading order logarithmic correc- 
tions to both the normal reflection amplitude and the 
AR amplitude which, in turn, give corrections to the con- 
ductance via the Landauer-Buttiker formula. In Sec. lIV| 
we obtain the RG equation for the NS junction and repro- 
duce the known fixed points using our approach. Then 
we derive the RG equation for the symmetric NSN junc- 
tion and obtain the RG flow between various fixed points 
and analyse the results of the study. Unlike, the renor- 
malization group flow of the NS junction, which does not 
lead to any non-monotonicity, we show that the inclusion 
of the CAR and direct tunneling through the supercon- 
ductor gives rise to a non-monotonic conductance as a 
function of the temperature. In Sec. El we present our 
results for specific cases of NSN . In Sec. lVIi we study the 
three— wire— superconducting— junction and show the ex- 
istence of a fixed point which is analogous to the Griffith's 
fixed poin t^^i^^ in the three— wire— normal— junction 
case. Finally in Sec. IVIH we present our summary and 
discussions. 



II. WIRG VIS-A-VIS BOSONIZATION 

Transport through a quantum scatterer (for instance, a 
simple static barrier or a dynamical impurity like Kondo 
spin) in a 1-D interacting electron gas is qualitatively 
different from its higher dimensional counterpart^. This 
is because, in 1-D, due to electron-electron interactions, 
the Fermi-liquid ground state is destroyed and the elec- 
trons form a non-Fermi liquid ground state known as 
as Luttinger Liquid^. The low energy dynamics of the 
1-D system is governed mainly by coherent particle-hole 
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excitations around the left and the right Fermi points. It 
is natural to use bosonic fields to describe these low ly- 
ing excitations. This can be done by re-expressing the 
original fermions using boson coherent state represen- 
tation^i^i^ which is referred to as bosonization. But 
this approach only allows for a perturbative analysis for 
transport around the limiting cases of SBS and WBS for 
the quantum impurity problem. On the other hand, if 
we start with a very weakly interacting electron gas, it 
is possible to do a perturbative analysis in the electron- 
electron interaction around the free fermion Hamiltonian, 
but treating the strength of the quantum impurity ex- 
actly. This allows us to study transport through the 
impurity for any scattering strength. The strength of 
this approach lies in the fact that even in presence of 
electron-electron interaction, one can use single particle 
notions such as the transmission and reflection ampli- 
tudes in order to characterize the impurity. The idea 
is to calculate correction to transmission and reflection 
amplitude perturbatively in the interaction strength. Of 
course, since we are working in 1-D, the perturbative 
correction turns out to be logarithmically divergent. To 
obtain a finite result, one has to sum up all such diver- 
gent contributions to the transmission and reflection am- 
plitudes to all relevant orders at a given energy scale. 
This was first done by Matveev et al. in Refs. [2^ 
and 27 in the context of a single (scalar) scatterer for 
both spin-less and spin-full electrons using the "poor 
man's scaling" approach In the spin-less case, it was 
shown that the logarithmic correction to the bare trans- 
mission amplitude (to first order in interaction parame- 
ter parameterized by a) was ST = 2 aTg (1 — To) ln(fcc?) 
and the explicit RG equation for transmission probabil- 
ity was dT/dl = ~2aT{l — T) where k was the momen- 
tum of the fermion measured from kp, d was a short 
distance cut-off and a was the interaction parameter 
given by a = ai — a2 with ai = V{0)/2TThvF and 
a2 = V(2kp)/2TThvF. The RG equation upon integra- 
tion gave the transmission probability as, 



T(L) = 
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1 - To + Toe-2"'] 



1 



To + To ay 



(1) 



Here, I = — In(fcd) = ln(L/d) where L is the length scale. 
I can also be measured as a function of the temperature 
by introducing the thermal length, Lt ~ {hvF)/{kBT). 
To is the bare transmission at the short distance cut- 
off, d. It is easy to see from Eq. [T] that for very small 
values of To, To can be neglected in the denomina- 
tor of the expression for T(L) leading to a pure power 
law scaling consistent with the power law known from 
bosonization in the WBS limit. Similarly for the spin- 
full electrons, it was shown that the parameter a in the 
power law gets replaced by a new parameter, f3 given by 
P = {92- 2gi)/TThvF where g2 = 52 (fc) and gi = gi{k) 
are momentum dependent functions or "running coupling 
constants". The momentum dependence of /3 here is in 
sharp contrast to the momentum independent a in the 



spin-less case. At high momentum, (or equivalently, at 
the short distance cut-off scale), gi{d) — V{2kF) and 
g2{d) = ^(0). Because of the extra logarithmic depen- 
dence coming from scaling of the interaction parameter 
itself (see Eq. [3] and Eq. |4] below), the expression for 
transmission^LSlj no longer shows a pure power law scal- 
ing even for small values of To. Instead T(L) is now given 

by 



T(T) 



(2q2-Qi) 



l-To + To[l + 2ailn(^)]^(^) 



d \ (2Q2-Qi) 



(2) 

using the length scale dependence of gi{L) and 32 (T) 
given bj*^ 
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V(2kF) 
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(4) 

Note that in the absence of electron-electron inter- 
action induced back-scattering (i.e., when V{2kF) = 
0), there is no correction to the power law behav- 
ior. Hence, bosonization, which ignores electron-electron 
back-scattering always results in power law behaviour. 
But, when electron-electron interaction induced back- 
scattering is included, the sign of 172 — '2gi can change 
under RG flow, and hence, there can be a qualitative 
change in the behavior of the conductance. The conduc- 
tance actually develops a non-monotonic dependence on 
the temperature; it first grows and then drops to zero. 
But, except for this non-monotonic behavior of conduc- 
tance for the spin-full case, there is no new physics which 
is found by studying the full crossover from WBS to SBS. 
In conclusion, both bosonization and WIRG methods pre- 
dict that for the single scatterer problem there are only 
two fixed points - (a) the perfectly back-scattering (no 
transmission) case is the stable fixed point and (b) the 
no back-scattering (perfect transmission) case is the un- 
stable fixed point. There are no fixed points with inter- 
mediate transmission. 

It was first shown by Lai et al. in Ref. i28i, using 
the WIRG approach that even though there are only 
two fixed points for the two— wire— junction, surpris- 
ingly enough, the three— wire— junction has a host of 
fixed points, some of which are isolated fixed points 
while others are one parameter or multi parameter fam- 
ilies of fixed points. It was also shown to be true 
for more than three wires. From this point of view, 
the physics of a two— wire— junction is different from 
its three— wire— counterpart. The three— wire— junction 
was also studied using bosonization and conformal field 
theory methods^'^'^'^, which confirmed some of the 
fixed points found using WIRG. It also gave some ex- 
tra fixed points which were related to charge fraction- 
alisation at the junction, and which could not be seen 
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within the WIRG approach. The WIRG method was 
further extended to compUcated systems made out of 
junctions of QW which can host resonances and anti- 
resonances in Ref. [H- The scahng of the resonances 
and anti-resonances were studied for various geometries 
which included the ring and the stub geometry. In par- 
ticular, it was shown that for a multiple— wire— junction, 
the RG equations for the full S'-matrix characterizing the 
junction take a very convenient matrix form, 

^ = -{SF^S - F) (5) 

where S is the scattering matrix at the junction and 
_F is a diagonal matrix that depends on the interaction 
strengths and the reflection amplitude in each wire. The 
advantage of writing the RG equation this way is that it 
immediately facilitates the hunt for various fixed points. 
All one needs to do is to set the matrix on the LHS of 
Eq. [5] to zero. This will formally provide us with all 
the fixed points associated with a given iS-matrix. This 
approach was further extended in Refs. I40ll4ll to study 
the multiple— wire— junction with a dynamical scatterer, 
i.e. a (Kondo) spin degree of freedom. The coupled 
RG equations involving the Kondo couplings, Jij as well 
as the 5'-matrices were solved. For different starting 
scalar S'-matrices, the RG flows of the Kondo couplings 
was studied. The temperature dependence of the con- 
ductances was shown to have an interesting interplay of 
the Kondo power laws as well as the interaction depen- 
dent power laws. Finally, the WIRG method was also 
extended to the case of NS junctio n^^i^° . In the vicinity 
of the superconductor, it is well-known that the system 
is described by holes as well as electrons^. Hence the S'- 
matrix characterizing the junction not only includes the 
electron channel but also the hole channel. Naively, one 
might expect that in the presence of particle-hole symme- 
try, the only effect of including the hole channel would 
be to multiply the conductance by a factor of two (in 
analogy with inclusion of spin and imposing spin up-spin 
down symmetry). However, it was shown^^'^*^ that in the 
vicinity of a superconductor, the proximity induced scat- 
tering potential that exists between electron and holes, 
also gets renormalized by electron-electron interactions. 
When this scattering is also taken into account, the cor- 
rection to the scattering amplitude to first order in the 
interaction parameter depends on (2g2 — gi) instead of 
{§2 — 2gi). It is worth stressing that this particular 
linear combination of the interaction parameters (g^'s) 
is independent of the scaling as the logarithmic factors 
{I — In(fcd)) in Eqs. [3] and [4] cancel each other. Hence, 
there is no non-monotonic behavior of the conductance in 
this case. The WIRG predicted only two fixed points, the 
Andreev fixed point (perfect AR) which turns out to be 
an unstable one and the perfectly reflecting fixed point 
which is the stable fixed point. The NS junction has also 
been studied using bosonization^. It is easy to check 
that the power laws resulting from bosonization agree 
with those obtained from the WIRG , when the electron- 
electron interaction induced back-scattering ( which is 



dropped in the bosonization method) is ignored. 

In this article, we apply the WIRG method to the su- 
perconducting junction of multiple quantum wires. We 
note that we now have two complications - (a) multiple 
wires are connected to the junction and (b) we have both 
electron and hole channels connected to the junction. So 
in this case, even for the NS junction, we have two spin 
channels as well as the electron and hole channels, so the 
scattering matrix is four component. For N wires, the 
scattering matrix is 4A^ x 4A^-dimensional. Although, we 
expect our method to work well even in this case, there is 
one caveat we must keep in mind. We have incorporated 
the effect of the superconductor as a boundary condition 
on the QW and neglected any internal dynamics of the 
superconductor itself. This should work reasonably well 
as long as we are studying transport at energies much 
below the superconducting gap. Our main result here 
is that the conductance across the junction depends on 
both gi and 52 and not on a a special combination 2g2—gi 
(as in NS case) which does not get renormalized under 
RG flow. Hence, the cancellation of the logarithmic terms 
in the effective interaction parameter is specific to the 
NS case and is not true in general. For N > 2 wires at- 
tached to a superconductor, we expect a non-monotonic 
form of the conductance. We also expect to get a host 
of fixed points with intermediate transmission and reflec- 
tion, knowledge of which can be of direct relevance for 
application to device fabrication of such geometries. 

III. SUPERCONDUCTING JUNCTION WITH 
N QUANTUM WIRES 

Let us consider multiple (N) quantum wires meeting 
at a junction on which a superconducting material is de- 
posited as depicted in Fig. [TJ The wires are parameter- 
ized by coordinates Xi, with the superconducting junction 
assumed to be at = 0. We consider a situation where 
the effective width 'a' of the superconductor between two 
consecutive wires is of the order of the phase coherence 
length of the superconductor (size of the Cooper pair). 
For our purpose, it is safe to ignore the finiteness of 'a' 
and effectively treat the junction of QW as a single point 
in space with an appropriate boundary condition. We 
parameterize the junction by the following quantum me- 
chanical amplitudes. There are two kinds of reflection 
amplitudes: the normal reflection amplitude (rusiSi) and 
the AR amplitude (rAusiSi) on each wire. In addition, 
there are two kinds of transmission amplitudes between 
different wire: the co-tunnelling (CT) amplitude (ty si sj) 
and the CAR amplitude (tAijsis)- The indices Si,Sj re- 
fer to the spin of incoming and outgoing particles. As we 
consider a singlet superconductor at the junction, spin 
remains conserved in all the processes mentioned above. 
Thus, the boundary conditions are parametrized by a 
AN X 4A^ scattering matrix for N quantum wires con- 
nected to a superconducting junction. 

Let us now consider the various symmetries that can 
be imposed to simplify the 4A^ x AN matrix. We impose 
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particle- hole symmetry, i.e., we assume that the reflec- 
tion and transmissions are the same for particles (elec- 
trons) and holes. Further, in the absence of a magnetic 
field, spin symmetry is conserved which implies that the 
various transmission and reflection amplitudes for spin 
up-down electrons and holes are equal. (This symme- 
try breaks down in the presence of magnetic fields, or 
in the case of ferromagnetic wires). Also, since we as- 
sume that all the wires, connected to the superconductor, 
are indistinguishable, we can impose a wire index sym- 
metry. (This symmetry again can be broken if we take 
some ferromagnetic and some normal wires attached to 
the superconductor). On imposing these symmetries, the 
S'-matrix for a two-wire system is given by 



5x 



with 



r t TA tA 

t r tA TA 

TA tA r t 

tA TA t r 



(6) 



Here r stands for normal reflection of electron or hole 
in each wire, and rA represents AR from electron to 
hole or vice-versa in each wire, t represents the elas- 
tic CT amplitude (t = ti2 = t2i) while tA represents the 
CAR amplitude {tA = tAi2 = tA2i)- For the spin sym- 
metric case, there are two such matrices, one for spin 
up electrons and holes and one for spin down electrons 
and holes. Note that this is the relevant S-matrix at 
energy scales (temperature and applied voltage on the 
wires) T, eVi <C A, where A is the superconducting 
gap energy. The competition between CT and CAR has 
been analysed before^-^ and also different ways of sepa- 
rating the contributions experimentally have been con- 
sideredi^. However, the effect of electron-electron in- 
teractions within the wires has not been considered for 
the NSN case. It is worth emphasizing here that if such 
NSN junctions are made out of 1-D systems like carbon 
nanotubes, then the effect of electron-electron interac- 
tions can influence the low energy dynamics significantly. 

The Landauer— Buttiker conductance matrix for the 
NSN case can be written, in the regime where ks T, eVi <^ 
A, asi 



11 _ Ga + GcA + GcT GcA — GcT Vi 

12 GcA — GcT Ga + GcA + Gct V2 

(7) 

The conductances here are related to the elements of the 
S'-matrix: Ga oc I^aP, Gct oc |ip and Gca oc |tAp- Ga 
is the conductance due to the AR that occurs at a single 
NS junction, whereas Gct and Gca are the conductance 
due to the elastic CT and CAR processes respectively, 
both of which involve transmissions between two wires 
and give contributions with opposite signs to the sub- 
gap conductance between the two wires, Gca — Gct ■ 
The generalization of this to > 2 is straightforward, 
and some details are presented in Sec. IVII 



IV. WIRG STUDY OF JUNCTIONS 

We study the effects of inter-electron interactions 
on the S'-matrix using the renormalization group (RG) 
method introduced in Re f. l27l. and the generalizations 
to multiple wires in Refs. I28ll29l The basic idea of the 
method is as follows. The presence of back-scattering 
(reflection) induces Friedel oscillations in the density of 
non-interacting electrons. Within a mean field picture 
for the weakly interacting electron gas, the electron not 
only scatters off the potential barrier but also scatters off 
these density oscillations with an amplitude proportional 
to the interaction strength. Hence by calculating the to- 
tal reflection amplitude due to scattering from the scalar 
scatterer and from the Friedel oscillations created by the 
scatterer, we can include the effect of electron-electron 
interaction in calculating transport. This can now be 
generalized to the case where there is, besides non-zero 
reflection, also non-zero AR . 

To derive the RG equations in the presence of Andreev 
processes, we will follow a similar procedure to the one 
followed in Ref. i28|. The fermion fields on each wire can 
be written as. 



(8) 



where i is the wire index, s is the spin index which can 
be tji and I,0 stands for outgoing or incoming fields. 
Note that ^i{x){^o{x)) are slowly varying fields on the 
scale of kp^ and contain the annihilation operators as well 
as the slowly varying wave-functions. For a momentum 
in the vicinity oi kp, the incoming and outgoing fields 
(with the incoming field on the i*'' wire) can be Fourier 
expanded in a complete set of states and the electron field 
can be written as 



*^s(a;) = 



dk 



dk 



+ tAdkse^'^'*' + *"'^'' + t*AhlE-'^^ + ^''^^ 



(9) 

where hus is the electron destruction operator and dks is 
the hole destruction operator. Note that we have cho- 
sen to quantize the fermions in the basis of the space 
of solutions of the Dirac equation, in the presence of a 
potential which allows for normal as well as Andreev scat- 
tering. We have also allowed for both incident electrons 
and holes. We find that (dropping a constant background 
density), 
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(10) 
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where the two terms corresponds to the density for elec- 
trons and holes respectively. Here we have also used 
the fact that due to the proximity of the superconduc- 
tor, the amplitude to create (destroy) a spin s electron 
and destroy (create) a spin s hole is non-zero — i.e., the 



Boguhobov amphtudes ( d\_gbl^ ) = 1 = ( bksdk-s ) 



sides the normal amplitudes 



dLdks 



be- 
= 1. 



(This is of course true only close to the superconduc- 
tor. We have checked that this gives the same result as 
solving the Boguliobov~de Gennes equation as done in 
Ref. l20j ) Hence, besides the density, the expectation val- 
ues for the pair amplitudes ( '^Is'^ls ) and its complex 
conjugate ( ^'i^^^s ) are also non-zero and are given by 
(dropping the wire index) 



and ( V'orV'/i ; 



A'KX 

Attx 



(11) 



So, we see that the Boguliobov amplitudes fall off as 1/x 
just like the normal density amplitudes. 

We now allow for short-range density-density interac- 
tions between the fermions 



1 



dxdy \ ^ Ps] V{x - y) 

5 = Ta 



5 = T,i 



(12) 

to obtain the standard four-fermion interaction Hamilto- 
nian for spin-full fermions as 



'H\, 



dx 
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T*OT*/T*OT 



.92(*lT*OT*OT*/T+*li*Oi*Oi*/i 



(13) 



where gi and (?2 are the running coupling constants de- 
fined in Sec. M (Eq. Eland Eq. H]). 

Using the expectation values for the fermion oper- 
ators, the effective Hamiltonian can be derived using 
a Hartree— Fock (HF) decomposition of the interaction. 
The charge conserving H F decomposition can be derived 
using the expectation values in Eq. [TO] and leads to the 
interaction Hamiltonian (normal) of the following form 
on each half wire, 



n 



N 



-i{g2 - 2ffi) 
47r 



dx r 

X 



r 



(*LT*/T+*Oi*/i 



(14) 

(We have assumed spin-symmetry i.e. = = ^0 
This has been derived earlier=. Using the same method. 



but now also allowing for a charge non-conserving HF de- 
composition with the expectation values in Eq. 1111 we get 
the (Andreev) Hamiltonian 



-»(gi +g2) 
47r 



dx 



,(15) 



Note that although this appears to be charge non- 
conserving, charge conservation is taken care of by the 
2e charge that flows into the superconductor every time 
there is an Andreev process taking place. 

The amplitude to go from an incoming electron wave 
to an outgoing electron wave under e"*'^!"'* (for electrons 
with spin) was derived in Ref. [2^ and is given by 



In(fcd) 



(16) 



where a — {g2 ~ 2gi) /2'KhvF and d was a short distance 
cut-off. Analogously, the amplitude to go from an incom- 
ing electron eip wave to an outgoing hole hout wave under 

£> '' ' '■in 



is given by 



= —I 



dk' 
2^ 



hout,s',fc' )( hout, s', A:' I U{^t |ein,s,A; 



92) r A 



dx 



-2ik X 



\Kuus\k') (17) 



where s ^ s' . Hence, the amplitude for an incoming 
electron to be scattered to an outgoing hole is given by 



a TA 



In(fcd) 



(18) 



where a' — {gi + g2)/'2.'KhvF- Note also that a and a' are 
themselves momentum dependent, since the giS are mo- 
mentum dependent. The amplitude for an incoming elec- 
tron to go to an outgoing electron on the same wire is gov- 
erned by the interaction parameter a = (g2 — 2gi)/2TrhvF 
which has the possibility of chaging sign under RG evolu- 
tion, because of the relative sign between 171 and g2. On 
the other hand, a' = (52 -I- gi)/2nhvF can never change 
its sign. 



A. NS Junction 

The amplitudes in Egs.fTBlandfTSlare corrections to the 
reflections of electrons from Friedel oscillations and from 
the pair potential respectively. We can combine them 
with the S'-matrix at the junction to find the corrections 
to the amplitudes of the S'-matrix. For an NS junction, 
there is only one wire coupled to the superconductor and 
the S'-matrix is just 2x2 for each value of the spin and 
is given by 



S = 



r TA 
TA r 



(19) 
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Here r is the normal refelction amplitude and ta is the 
Andreev reflection amplitude. So we only need to com- 
pute the corrections to r and ta in this case. 

We find that there are five processes which contribute 
to the amplitude to first order in the interaction pa- 
rameter. These are illustrated in Fig. [2l 




FIG. 2: The processes that contribute to the amplitude for an 
incoming electron to transform to an outgoing hole. Note that 
all the processes shown here are to first order in the interaction 
parameters since they only involve a single scattering from a 
Friedel oscillation or the pair potential. Process (c) involves 
scattering from a pair potential before the electron reaches 
the junction. The remaining processes involve two reflections 
from the junction and a scattering from the Friedel oscillation 
or the pair potential. In the diagrams, ^ = ^ar*\n{kd), 
rj — — r\ \n{kd) and r)' — ta \vdkd). 

Adding all the contributions, we obtain the change in 
the AR amplitude ta that takes an incoming electron to 
an outgoing hole given by 



8r I 



a' 

y 



TA 



In(fcrf) 



-f a |rp TA \n{kd) 



(20) 



in agreement with Ref. |43|. For an incoming electron 
refiected back as an electron, we find the small correction 
in the amplitude 6r given by^li^^ 



Sr = —a'\rA\'^r\n{kd) 
+ ? [rlr* - r(l 



In(fcd) (21) 



We replace — In(fcd) by dl using the "poor man's scaling" 
approach^^ to obtain the RG equation for ta as 

^ = -^[rA^r*A{r' +r\)]-a\r\\A (22) 

Using the unitarity of the S"- matrix (| '"a P + | ^"^ | = 1 
and r^r + tat* — 0), we can simplify the RHS of the 
above equation to obtain 



^ = -{a + a')rA{l-\rA\' 



(23) 



Note that the combination a + a' ~ (2 .92 — .9i)/27r?iUi? 
which appears in the RG equation does not flow under 
RG. This can be seen from Eqs. [3] and [4] which shows 
that {2g2-gi)/27rhvF = {2V{0)-V{2kF))/2T:hvF- This 
means that r and either monotonically increase or 
decrease as a power law depending on the sign oi a + a' . 
From Eq. [531 we also observe that |r^| = and |r^| = 1 
correspond to the insulating and the Andreev fixed points 
of the NS junction respectively. One can easily see from 
the RG equations that |r^| = is a stable fixed point and 
jr^l = 1 is an unstable fixed point. 



B. NSN Junction 

In this subsection, we shall consider an NSN junction. 
Here in addition to the two reflection channels, we also 
have two channels for transmission - the direct transmis- 
sion of an electron to an electron through CT process and 
the transmission of an electron to a hole via CAR. These 
processes are depicted in Fig. [3] The ^-matrix at the 
junction is 8 x 8 in this case as given in Eq. [6l The num- 
ber of processes that contribute in this case is thirty four, 
since we also need to include terms that transmit elec- 
trons or holes through the junction. For instance, for the 
renormalization of the AR term, besides the terms cor- 
responding to the NS junction, we also have to include 
processes in which the electron is incident from wire 1, 
goes through the junction to wire 2, Andreev reflects from 
the pair potential on wire 2 and then comes back through 
the junction, as shown pictorially in Fig. SJc). 

Collecting all the nine processes that contribute to first 
order in a and a' to the reflection amplitude, we find that 



dr 
'dl 



a'{r\rAf + r*AtAt) (24) 








N 


S 


N 

> 



FIG. 3: Electron CT with bare amplitude t is shown in the 
left plot and CAR with bare amplitude Ia is shown in the 
right plot. 
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Similarly, adding up the contributions from the nine pro- 
cesses that contribute to r^, we find that 



drA 
dl 



a(|rprA + ttAr*) 



+ - {r^ + r\ + e + t\)r\) 



(25) 
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FIG. 4: The extra processes that contribute to the amplitude 
for an incoming electron to transform to an outgoing hole on 
the same wire, due to the second wire. Processes (a) and (b) 
are transmitted to the second wire and reflected by the Friedel 
oscillation whereas (c) and (d) are transmitted to the second 
wire and reflected by the pair potential. 

Moreover, here besides the reflection parameters, we 
also need to compute the renormalizations of the trans- 
missions to first order in a and a' . The RG equations for 
t and tA are also obtained by considering all possible pro- 
cesses that ultimately have one incoming electron and one 
outgoing electron (for t) and one incoming electron and 
one outgoing hole (for Ia) and are either reflected once 
from the Friedel potential or the pair potential. They are 
found to be 



dt 
Jl 



dtA 
dl 



adrpt + r^TAtA) 
a'{\rA\^t -f rr\tA) 
- a{r*rAt + |rp i^) 
a'{rtr\ + |r^pt^) 



(26) 



(27) 



Just as was done for the normal junction (EqlS]), we 
can express the RG equations for the superconducting 
junction in a compact matrix form^^. 



where the matrix S is given in Eq[6] and F depends on 
the interaction parameters a ~ {§2 — 2gi) /2TThvp and 
a' — {gi + g2)/2'Khvp. F is non-diagonal matrix (unlike 



the case in Ref. '2g) and is given by 



ar 
2 





— a rA 
2 








ar 
2 





-a'rA 
2 


a'rA 
2 





ar 
2 








-a'rA 
2 





cxr 
2 



(29) 



^=F-SF^S 
dl 



(28) 



It is easy to check that all the RG equations are repro- 
duced from the matrix equation. The matrix form also 
makes the generalization to N wires case notationally 
simple and makes the search for various fixed point much 
easier. This will be discussed in the last section. But note 
that these equations have to be augmented by Eqs.[3]and 
|4]to get the full set of RG equations. 

Let us now look at some of the fixed points of the 5*- 
matrix. Clearly, the fixed points occur when F — SF'^S = 
or when FS'^ is hermitian. There are several possibili- 
ties and we list below some of them. 

Case I: Any one of the four parameters is non-zero 

(a) t = \, r = TA = tA = fully transmitting fixed 
point (TFP) 

(b) r = 1, rA = tA = t = Q fully reflecting flxed 
point (RFP) 

(c) rA — 1, r — t — tA — 0, fully Andreev reflecting 
fixed point (AFP) 

(d) tA = r = t = TA = 0, fully crossed Andreev 
reflecting fixed point. (CAFP) 

Case II: Any two are non-zero 
When both r and rA are zero, the RHS of the 
RG equations identically vanishes as both the 
Friedel oscillation amplitude as well as the pair po- 
tential amplitude in the wire become zero. Hence 
any value of t and tA remains unrenormalized under 
RG. 

Case III: Any three of them are non-zero 
We did not find any fixed point of this type. 

Case IV: All four of them are non-zero 
Here, we get a fixed point when ri = r2 = t = tA = 
1/2 and r^i — M2 = ~l/2. This is the most sym- 
metric S'-matrix possible for the NSN case. Since 
it is a symmetry-dictated fixed point with interme- 
diate transmission and reflection, we shall refer to 
it as symmetric fixed point (SFP). 

We will study the RG flows near some of these fixed points 
in the next section. 



C. FS, FSF and FSIM Junctions 

We can also consider junctions where one or more of 
the wires are spin-polarised, with Fermi distributions for 
the spin up and down electrons being different. As long 
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as at least one of the wires is ferromagnetic, the spin up- 
spin down symmetry of the system is broken. This means 
that we can no longer impose S-^ = on the ^-matrix 
parametrising the scattering as we had in Eq|Bl We now 
need to choose an S'-matrix with indices t and | denot- 
ing the spin. For the FSN case (and the FSF case where 
the ferromagnets on the two sides are not identically po- 
larized) the wire index symmetry is also broken. Hence, 
the S'-matrix chosen must also break the wire-index sym- 
metry. Note that for the ferromagnetic wire, the ampli- 
tude to destroy a spin s electron and create a spin s hole 
cannot be non-zero, even in the proximity of the super- 



conductor. The Boguhobov amplitudes ( dli._^bli^^ ) and 
( biksdik-s ) decay exponentially fast (with a length scale 
set by the ferro— anti-ferro gap) in the ferromagnetic wire. 
So, in our S-matrix, ta is zero and there is no pair po- 
tential due to the proximity effect in ferromagnetic wire. 
Also as mentioned earlier, we must keep in mind that 
the influence of the ferromagnet on the spectrum of the 
superconductor has to be negligibly small. This will be 
true only if the superconductor is large enough. Hence, 
for such junctions, the rcnormalization of the S'-matrix 
is only due to the Friedel oscillations. Also note that in 
these wires, since the bulk does not have both the spin 
species, gi and 52 do not get renormalized. All the cases 
mentioned above will therefore involve the full 4N x AN 
S-matrix since there is no reduction in number of inde- 
pendent elements of the S-matrix which can occur when 
symmetries are imposed. 



D. Three Wire— Junction A Beam Splitter 



ing in transport. Then the S-matrix is given by 



S: 












rA 


tA 


t'A 











tA 


TA 


t'A 











t'A 


t'A 


r'A 


rA 


tA 


t'A 











tA 


TA 


t'A 











t'A 




r'A 












(30) 



Here, ta and Ia and are real parameters which sat- 
isfy^^ 



tA 



= 1 + , 
= -l-'ivA 



and 



= V(-2M)(l + r^) 

- 1 < < 



(31) 



by unitarity. Using Eq. the simplified RG equation 
for the single parameter ta is given by 



dVA 

dl 



a' [rA{l + rA){l + irA)] (32) 
So, within the real parametrization we have two unsta- 





r =-1/3 
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— * --^-^ 

r =0 
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In this subsection, we consider the standard beam 
splitter geometry comprising of a superconductor at the 
junction of three quantum wires. In this case, we show 
that there is a fixed point that is analogous to the An- 
dreev fixed point of the NS junction. The S-matrix rep- 
resenting this fixed point is symmetric under all possi- 
ble permutations of the three wires and allows for the 
maximum Andreev transmission (in all channels simul- 
taneously within unitarity constraints). The S-matrix 
is given by = —1/3 and = = 2/3 with 
r = t = t' ^ Q. We refer to this fixed point as the 
Andreev- Griffith's fixed point (AGFP)il. 

For an analytic treatment of this case, we will consider 
a simplified situation where there is a complete symmetry 
between two of the wires, say 1 and 2, and the S-matrix 
is real. In addition, the elements of the S-matrix corre- 
sponding to transmission or reflection of an incident elec- 
tron (hole) to a reflected or transmitted electron (hole) 
are set to zero so that only Andreev processes participat- 



FIG. 5; Schematic representation of the situation where a 
three— wire— junction is hooked to the stable fixed point, 
AGP P. An incident electron in one wire is either reflected back 
as a hole in the same wire or is transmitted as a hole in an- 
other wire along with the addition of the two electrons into 
the superconductor forming a Cooper pair. The direction of 
RG flow from two unstable fixed points to the stable fixed 
point (AGFP) is also depicted on the bottom left side of the 
diagram. 

ble fixed points, given by r^i = and r^i = — 1 and a 
stable fixed point given by = — 1/3. The = fixed 
point corresponds to a situation where there is perfect 
CAR between wires 1 and 2 and wire 3 gets cut off from 
the remaining two wires (labelled by 1 and 2) and is in the 
perfect AR limit with the superconductor. The = — 1 
fixed point corresponds to a situation where all the three 
wires are disconnected from each other and are in perfect 
AR limit individually with the superconductor. The third 
fixed point given by = —1/3 corresponds to a perfect 
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Andreev limit of the three wire junction where an in- 
cident electron is either Andreev-reflected into the same 
wire as a hole or is transmitted as a hole via CAR into an- 
other wire. This is essentially the AGFP. It is very inter- 
esting to note that the original Griffith's fixed point was a 
repulsive fixed poin t^^i^^ whereas the AGFP is an attrac- 
tive fixed point. This can be understood as follows. Here, 
there is no scattering from the Friedel oscillations as the 
junction is assumed to be reflection-less, whereas there 
exists a proximity induced pair potential, which induces 
an effective attractive interaction between the electrons. 
Hence, the physics is very similar to the well-known Lut- 
tinger Liquid physics, which says that for attractive in- 
teraction between the electrons, back-scattering is an ir- 
relevant operator. Hence the stable fixed point here will 
be the one which will have maximal transmission between 
the wires. So, it is not surprising that the AGFP turns 
out to be a stable fixed point. Thus, for a rcflcction-lcss 
junction, we have found a stable fixed point with inter- 
mediate transmission and reflection. 




10 20 30 40 



L 



FIG. 6: Conductance of the NS junction is plotted in units 
of e^//i as a function of the dimensionless parameter / where 
I = ln{L/d) and L is either Lt = hvp/ksT at zero bias or 
Lv ~ hvp/eV at zero temperature and d is the short distance 
cut-off for the RG flow. The three curves correspond to three 
different values of V'(O) and V{2kF). 



V. RESULTS 

In this section, we will consider various physical cases 
and see what the RG flows mean for the conductances in 
each case. 



A. NS Junction 

First, we give the results for the NS junction, just to 
contrast with the results of the NSN junction. Here, we 
have only two parameters, r and va- The conductance 
occurs only due to the AR amplitude, va which obeys 
the RG equation given by Eq. [201 As mentioned already, 
there is no flow of the particular linear combination of 
the interaction parameters 2g2 — gi that occurs in the 
equation and the RG flow of the conductance is there- 
fore monotonic. The conductance as a function of the 
length scale for different interaction parameters V{Q) and 
V{2kp) is plotted in Fig. [SI L here simply denotes the 
length at which the RG is cut-off. So if we take very long 
wires ^ Lt, then the cut-off is set by the temper- 
ature, and the plot shows the variation of the conduc- 
tance as a function of Lt starting from the high tem- 
perature limit, which here is the superconducting gap A. 
We observe that as we lower the temperature, the An- 
dreev conductance decreases monotonically. Also it was 
established in Refl20l that the power law scaling of con- 
ductance ( It'aP) calculated from WIRG and bosonization 
were found to be in agreement with each other for the 
limiting cases of |rAp = 1 and |rAp = ( which are the 
only limits where bosonization results are valid) provided 
effects due to electron-electron induced back-scattering in 
the wires is neglected. 



B. Ballistic NSN Junction 

In this subsection, we consider the case of a reflection- 
less ballistic junction between the superconductor and 
the two wires, i.e. r = 0. This implies that the renormal- 
ization of the S'-matrix due to the Friedel oscillations is 
absent. The only renormalization is due to reflections 
from the proximity effect induced pair potential. Let us 
now consider various interesting cases: 

(a) r = 0, = 0, t 7^ 0, <yi 7^ 0. In this case, since 
we have both r — Q and r^i = 0, there is no RG flow 
of the transmission and the conductance is frozen 
at the value that it had for the bare 5'-matrix. The 
most interesting situation in this case arises when 
t = tA- For this case, the probability for an incident 
electron in one wire, to transmit in the other wire 
as an electron due to t or as hole due to tA is equal, 
leading to perfect cancellation of charge current. 

(h) r = 0,t = 0,rA ^ 0,tA ^ 0- For this case, one 
can easily check from the RG equations (Eqs. [Ml ^ 
\T7]\ that if we start our RG flow with the given 
parameters at high energies, then the value of r, t 
remain stuck to the value zero under the RG flow. 
Hence, in this case the two parameter subspace 
rAy^O,tAy^O remains secluded under the RG flow. 
The RG equation for tA is given by 

^ = a'tA{l~\tAf) (33) 

The above equation can be integrated to obtain an 
expression for CAR probability {Ta — I^Ap), 
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rpO 


[1 


+ 2ailn(^)]i(^)-(2"2-ai)' 






[l + 2ailn(^)]i(|)-(2"2-«i) 



T4 and R\ are the CAR and AR probabilities at 
the short distance cut-ofF, L = d. We notice that 
the RG equation and its solution are very similar to 
that for the single scatterer problem^ apart from 
a sign difference on the RHS of the equation and 
the dependance of the interaction parameter a' on 
gi and 52 ■ This implies that even if we start with a 
small crossed Andreev transmission across the junc- 
tion, the RG flow will take us towards the limit of 
perfect transmission. This is in sharp contrast to 
the normal transmission across a single scatterer. 
For the single barrier problem, the equation for the 
RG flow of t was by 

I = -at(l-|t|2). (35) 

Hence, t = was the stable fixed point. But if 
the electron-electron interactions had been attrac- 
tive, then the sign on the RHS would have been 
positive and t — 1 would have been the stable 
fixed point. Thus, the RG flow of Ia for the case 
when rA 7^ 0,t = r = 0, and repulsive interac- 
tions, is very similar to the RG flow for t when 
r 7^ 0, = = but with attractive interac- 
tions. In both cases transmission is relevant and 
t = 1 and Ia = ^ are the stable fixed points. 
On the other hand the RG flow of Ia for the case 
offA 7^ 0,t = 0,r = and attractive electron- 
electron interaction (t^(0), V{2kp) < 0) in the wire 
is very similar to the RG flow for t for the case 
r 7^ 0,rA — Q,tA — and repulsive electron- 
electron interaction {V{Q),V{2kp) > 0). In both 
cases transmission is irrelevant and t — Q and 
tA — Q are the stable fixed points. At an intu- 
itive level, one can perhaps say that even if we 
start with repulsive inter-electron interactions, the 
proximity-induced pair potential leads to a net at- 
tractive interaction between the electrons, which is 
responsible for the counter-intuitive RG flow. 

Also notice that while solving the above RG equa- 
tion for Ia, we have to take into account the RG flow 
of the interaction parameter (a') itself. This will 
lead to non-power law (non Luttinger) behavior for 
the conductance close to ji^j ~ or ji^j ~ 1. 
It is worth pointing out that the non-power law 
part appearing in Eq. [M] is identical to Ref. H^, 
even though the interaction parameter for their 
case was proportional to (72 — 2(7i and for our case 
it is (72 + 31 ■ But of course this will not lead to any 
non-monotonic behavior as a' can not change sign 
under RG flow. So the stable fixed point for this 
case is the CAFP. 



(c) r = 0, = 0, 7^ 0, t 7^ 0. This case is iden- 
tical to the case (b) discussed above except for the 
fact that we have to replace Ia in the previous case 
by t. In this case also the two parameter subspace 

7^ 0, i 7^ remains secluded under RG flow. The 
RG equation for t is given by 

I = a't{l~W) (36) 

Here also, t = 1 remains the stable fixed point and 
i = is the unstable fixed point. 

(d) r = Q,t ^ Q,rA Q,tA ^ 0. In this case if 
we start from a symmetric situation, i.e. t = Ia, 
we can see from the RG equations in Eqs. [26] and 
[27] that both t and tA have identical RG flows. So, 
the sub-gap conductance G = Gca — Gct van- 
ishes identically and remains zero through out the 
RG flow. Hence this ^-matrix can facilitate pro- 
duction of pure SC— if we inject spin polarized elec- 
trons from one of the leads as the charge current 
gets completely filtered out at the junction. 



C. Ballistic FSF Junction 

Here, we consider the case where both the wires are 
spin polarized. In this case we have two interesting 
possibilities, i.e. either both the wires have aligned spin 
polarization or they have them anti-aligned. In either 
case the Andreev reflection amplitude is zero on each 
wire due to reasons explained earlier. 



(a) When the two wires have their spins aligned, 
t 7^ 0, but Ia = because for CAR to happen we 
need up (down) spin polarization in one wire and 
down(up) spin polarization on the other wire which 
is not possible in this case. 

(h) When the two wires have their spins anti- 
aligned, i = 0, but Ia ^ ^ because the up(down) 
electron from one wire can not tunnel without flip- 
ping its spin into the other wire. As there is no 
mechanism for flipping the spin of the electron at 
the junction, such processes are not allowed. 

Hence these two cases can help in separating out and 
measuring amplitudes of the direct tunneling process and 
the CAR process experimentally-'^^. Both these are exam- 
ples of case II, since they have both r — and = 0. 
In this case, neither t nor Ia change under RG flow and 
hence conductance is not influenced by electron-electron 
interaction at all. 
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FIG. 7: Conductance Gca of the NSN junction is plotted 
(when the two leads have anti-parallel spins) in units of 
/h as a function of the dimensionless parameter I where 
I = ln{L/d) and L is either Lt = hvp/ksT at zero bias or 
Lv = hvp/eV at zero temperature and d is the short dis- 
tance cut-off for the RG flow. The three curves correspond to 
three different values of V(0) and V{2kF)- The inset shows 
the behavior of the same conductance for fixed values of a. 



D. Non-ballistic NSN Junction 

1. Without AR on individual wires 

Here we consider an NSN junction with finite reflection 
in each wire and no AR in the individual wires. So the 
renormalization of the S-matrix is purely due to Friedel 
oscillations and there are no contributions coming from 
scattering due to the pair potential. Below we discuss 
two cases : 



(a) rA — 0,t ~ Q,r ^ Q,tA 7^ 0. This is an exam- 
ple of case II mentioned in Subsection. IIVBI The 
RG equations (Eqs. [^ill?f|) predict that rA,t will 
remain zero under the RG flow and r, Ia form a se- 
cluded sub-space. The RG equation for this case is 
given by 



To contrast, we also show in the inset, the behavior 
when the renormalization of a in not taken into ac- 
count. Thus, it is apparent from the plot that the 
non-monotonicity is coming solely from the RG evo- 
lution of a. The inset and the main graph, both 
start from the same value of tA- Even though this 
case is theoretically interesting to explore, its ex- 
perimental realization may not be viable. This is 
because of the following reasons. Here we have 
r^i = on both wires, which can only happen if the 
wires are ferromagnetic. However, we also know 
that if the wires are ferromagnetic, there is no scal- 
ing of a parameter and hence there will be no inter- 
esting non-monotonic trend in the conductance. So 
it is hard to flnd a physical situation where = 
and at the same time, there is renormalization of 
the interaction parameter a. Lastly note that the 
conductance Gca is negative. The process respon- 
sible for the conductance, (i.e. CAR), converts an 
incoming electron to an outgoing hole or vice- versa, 
resulting in the negative sign. 



(b) rA — 0,tA = 0,r ^ 0,t ^ 0. This case is iden- 
tical to the previous case with the replacement of 
tA by t. Fig. [8] shows the the CT conductance Gct 
as a function of the length scale. It shows a simi- 
lar non-monotonic behavior with positive values for 
the conductance. The inset shows the behavior of 
Gct when the renormalization of a in not taken 
into account. 



V(0)=0.1=V(2kp) 
V(0)=0.2=V(2kp) 
V(0)=0.3=V(2k ) 
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^ = -atAil-\tA\') 



(37) 




Note the change in sign on the RHS with respect 
to the RG equation for tA (Eq. [33]) for the ballis- 
tic case. This change in sign represents the fact 
that the ballistic case effectively represents a situ- 
ation corresponding to attractive electron-electron 
interaction while this case corresponds to a purely 
repulsive electron-electron interaction. 

Fig. [7] shows the behavior of conductance {Gca) 
for this case. The conductance in the main graph 
shows a non-monotonic behavior. 



FIG. 8: Conductance Gct of the NSN junction (when the 
two leads have parallel spins) in units of /h as a function 
of the dimensionless parameter I where I — ln{L/d) and L 
is either Lt = hvp/ksT at zero bias or Lv = hvp/eV at 
zero temperature and d is the short distance cut-off for the 
RG flow. The three curves correspond to three different values 
of V'(O) and V{2kF)- The inset shows the behavior of the same 
conductance for fixed values of a. 
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2. With AR on individual wires 

(a) rA 0,tA 7^ 0,r 0,t ^ 0. This is the most 
interesting case, where both r and va are non-zero, 
and we get an interplay of the effects due to scatter- 
ing from Friedel oscillations and from proximity in- 
duced pair potential. Here, all the four parameters 
are non-zero and flow under RG, as do the interac- 
tion parameters a and a'. An example where the 
system starts in the vicinity of the unstable fixed 
point SFP (as mentioned in Case IV in the Sub- 
section. IIVBP is shown in Fig. [31 The NSN con- 
ductance here is defined as Gnsn — Gca — Gct- 
Here also we observe a strong non-monotonicity in 
the conductance which comes about due to inter- 
play of the electron and the hole channels, which 
contribute to the conductance with opposite signs, 
coupled with the effects from the RG flow of the 
interaction parameters. 



E. Non-ballistic FSN Junction 

In this case, for the ferromagnetic wire = 0, but 
for the normal wire ta has a finite value. As explained 
earlier, the interaction parameters a and a' on the fer- 
romagnetic side do not renormalize, whereas they do on 
the normal side. Hence, even if we start from a situation 
where the interaction parameter a and a' are symmet- 
ric for the two wires, RG flow will always give rise to 
an asymmetry in the interaction strength. Therefore, it 
becomes a very interesting case to study theoretically. 
The S'-matrix for this case has neither spin up-spin down 
symmetry, nor the wire index (left-right for two wires) 
symmetry. Only the particle-hole symmetry can be re- 
tained while parameterizing the S'-matrix. This case gets 
very complicated to study theoretically because the mini- 
mum number of independent complex-valued parameters 
that are required to parameterize the S'-matrix is nine 
as opposed to four in the NSN case. These are given by 

„11 „22 22 fl2 f21 22 22 j.12 j j.21 

'TT' TT' il' ^TT' Ai[^ 'ATT' ^ii' TT' '^"'^^tT' ^^^^^ 
1(2) is the wire index for the ferromagnetic (normal) wire 
while, t and [ are the respective spin polarization indices 
for the electron. 

So, the minimal S-matrix representing the FSN junc- 
tion is given by 



r 
t' 



t'A 





t tA 
r' TA 
r" t'^ 



tA r 

r'A 
VA t' 



t 



r" 



(38) 



The RG equations for the nine independent parameters 
are given in Appendix |^ We write down a represen- 
tative S'-matrix which satisfies all the constraints of the 
FSN junction and unitarity, and study its RG flow nu- 
merically by solving the nine coupled differential equa- 
tions. The modulus of the S-matrix elements are given 



by 



„ii| 



ttI 



I^attI 



nil 



•22 



TTl 



= t 



TTl 



= 1/V3 and the correspond- 
ing phases associated with each of these amplitudes are 
7r/3, TT, 0, — 7r/3, 0, 7r/3, 0, tt, — 7r/3 respectively. Here also 
we observe a non-monotonic behavior of conductance, 
Gfsn sls a function of I as shown in Fig. (5] 



F. Non-ballistic FSF junction 



Here, we will consider the case where both the wires 
are spin polarized. This case is similar to the ballistic 
case. Since here r 7^ and ta = 0, we will have the usual 
Friedel oscillations and the conductance will go to zero as 
a power law. Here again we have two instructive cases: 



(a) When the two wires connected to the supercon- 
ductor have their spin polarization aligned, i.e. t ^ 
0, but tA = and 



(b) When the two wires have their spin polarization 
anti-aligned, i.e., t = but tA 7^ 0. 



Both these are examples of case II of Subsection. IIVBI 
Either t or tA need to be zero in the two cases mentioned 
above. Hence the parameters which are zero will remain 
zero under RG , while the non-zero parameters will 
flow according to Eqs. [26] and [27] respectively. The 
conductances are the same as in the NSN case except 
that the interaction parameters cannot flow now. This 
has already been plotted in the insets in Figs. [7] and [5] 
Since the electrons are now effectively spin-less, a and 
a' do not flow, and we get a monotonic fall-off of the 
conductance in both the cases. 

The results of this section are summarized in the table 
below. 
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L L 

FIG. 9: Left: Conductance of the NSN junction Givsiv ~ [iAp~|ip is plotted in units of 2e^//i as a function of the dimensionless 
parameter I where I = ln{L/d) and L is either Lt = hvp/kBT at zero bias or Ly = hvp/eV at zero temperature and d is the 
short distance cut-off for the RG flow. The three curves correspond to three different values of V{Q) and V{2kF)- Right: Charge 
conductance Gfsn = \tA\^ — \t\^ is plotted for FSN case as a function of the dimensionless parameter / where / — ln{L/d) and 
L is either Lt = hvp/ksT at zero bias or Lv ~ hvp/eV at zero temperature in units oi e^/h and d is the short distance cut-off 
for the RG flow. The three curves correspond to three different values of V(0) and V{2kF)- 



t 


tA 




r 


Stability 


Intermediate fixed point 


Relevant physics 











1 


Stable 


X 


RFP 








1 





Unstable 


X 


AFP 





1 








Unstable 


X 


CAFP 


1 











Unstable 


X 


TFP 


1/2 


1/2 


-1/2 


1/2 


Unstable 




SFP, Non-monotonic charge current 


e^'^i sin 6* 


e'-^^ cos e 








Marginal 




Pure spin current when t — Ia 



VI. GENERALIZATION TO THE CASE OF 
THREE WIRES 

In this section, we consider the case of three wires con- 
nected to a superconductor. We assume that all the wires 
are connected within the phase coherence length of the 
superconductor. Hence, CAR can occur by pairing the 
incident electron with an electron from any of the other 
wires and emitting a hole in that wire. The conductance 
matrix can hence be extended for three wires as 



(39) 



h 




Gr 11 Gt 12 Gt 13 




Vi 


h 




Gt 12 Gr 22 Gt 23 




V2 


h 




Gt 13 Gt 23 Gr 33 




V3 



with Gr 



G, 



" J2j{GcAij + GcTij) and Gtij = 
GcAij — GcTij and the generalization to N wires is 
obvious. Note that the conductances GcAij — GcAji 
and GcTij — Gct 



The relations of the conduc- 
tances to the reflections and transmissions is obvious, for 
e.g., Gau oc Ir^iiP as before while Gcnj oc \tij\'^ and 
GcAij oc The RG equations for the three wire 

case can be written using the matrix equation as given 
in Eq. [6] except that the S matrix is now 12 x 12 di- 
mensional. For a system with particle-hole, spin up-spin 



(40) 



down and wire index symmetry, the S matrix is given by, 

r t t' ta tA t'^ 
t r t tA TA tA 

t' t r t'^ tA Ta 

TA tA t'^ r t t' 

tA rA tA t r t 

t'A tA TA t' t r 

where we have chosen six independent parameters, with 
ti2 = t2i = t23 = t^2 = t and ii3 = t^i = t' and similarly 
for the CAR parameter tA- The F matrix now generalizes 
to 



F 



OLT 
2 








-a TA 
2 











ar 
2 








-a'rA 
2 











ar 
2 








-a'rA 
2 


— a'rA 
2 








ar 
2 











-a'vA 
2 








ar 
2 











2 








ar 
2 



(41) 

and the RG equations for the six independent parame- 
ters are given in Appendix [B] There exists possibility 
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of many more non-trivial fixed points in this case. For 
instance, the AGFP, as mentioned in Subsection. IIVDI 
As discussed in Subsection. IIVD| for the reflection-less 
case with symmetry between just two wires, this compli- 
cated S'-matrix takes a very simple form, which can be 
dealt with analytically. Within the sub-space considered 
we found that the AGFP was a stable fixed point. In 
Fig. [ini we show the RG flow of from two different 
unstable flxed points to the stable AGFP. 

The possibility of experimental detection of such a non- 
trivial flxed point with intermediate transmission and re- 
flection is quite interesting. From this point of view, the 
AGFP is a very well-suited candidate as opposed to its 
counterpart, the Griffith's fixed poin t^^'^^ . For a normal 
junction of three 1-D QW, the S'-matrix corresponding to 
r = —1/3, t = 2/3 is a fixed point (Griffith's fixed point), 
where r and t are the reflection and the transmission for a 
completely symmetric three wire junction. Even though 
it is an interesting flxed point, it turns out to be a re- 
pulsive one and hence the possibility of its experimental 
detection is very low. On the contrary, the AGFP, being 
an attractive fixed point, has a better possibility of be 
experimentally measured. The main point here is that 
even if we begin with an asymmetric junction, which is 
natural in a realistic experimental situation, the effect of 
interaction correlations are such that as we go down in 
temperature, the system will flow towards the symmetric 
junction. This can be inferred from the results shown in 
Fig.Iini 
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FIG. 10: |tAp is plotted as a function of dimensionless pa- 
rameter I where I — ln{L/d) and L is either Lt ~ hvp/ksT 
at zero bias or Lv = hvp/eV at zero temperature and d is 
the short distance cut-ofF for the RG flow. The three curves 
correspond to three different values of V{0) and V{2kF)- The 
set of curves in the top represent the RG flow of \tA P when the 
starting point is in the vicinity of ta = fixed point while, 
the set of curves in the bottom half represent the RG flow of 
|tA|^ with starting point close to ta = — 1. 



VII. SUMMARY AND DISCUSSIONS 

To summarize, in this article we have studied trans- 
port through a superconducting junction of multiple 
1-D interacting quantum wires in the spirit of the 
Landauer— Buttiker formulation. Using the WIRG ap- 
proach we derived the RG equations for the effective S'- 
matrix and obtained the various flxed point S-matrices 
representing the junction. In contrast to earlier RG stud- 
ies, here, we had to include both particle and hole chan- 
nels due to the proximity effect of the superconductor. 
Our study led to the flnding of a novel flxed point with 
intermediate Andreev transmission and reflection even 
in the case of NSN junction (SFP). However, it turns 
out to be unstable flxed point and hence experimentally 
inaccessible. We found that transport across the super- 
conducting junction depends on two independent inter- 
action parameters, (a = (52 — '^gi) I'^t^^vf) which is due 
to the usual correlations coming from Friedel oscillations 
for spin-full electrons and (a' = (52 + ffi)/27rSwi?), which 
arises due to the scattering of electron into hole by the 
proximity induced pair potential in the QW. We com- 
puted the length scale (or temperature) dependance of 
the conductance taking into account interaction induced 
forward and back-scattering processes. We found a non- 
monotonic dependence of the conductance for the two- 
wire NSN superconducting junction in contract to the 
NS junction where the dependence is purely monotonic. 

When more than two wires are attached to the super- 
conductor, we found even more exotic flxed points like 
the AGFP which happens to be a stable flxed point for 
a reflection-less symmetric junction thus increasing its 
chances of being experimentally seen. But reflection is a 
relevant perturbation and asymmetry between different 
wires is likely to be a relevant parameter—. Hence, the 
RG flow due to these perturbations will take the system 
to the RFPin the low energy limit. So, the experimen- 
tal detection of AGFP flxed point will critically depend 
on how efflciently the conditions of reflectionlessness and 
symmetry can be maintained in the experimental setup. 
If the reflection and asymmetry are reduced to a large 
extent, then as we cool the system, the S-matrix at the 
j miction is expected to flow very close to the AGFP flxed 
point. But ultimately, the RG flow will lead to enhance- 
ment of any initial small value of reflection and asym- 
metry and we will finally flow to the disconnected fixed 
point in the T = limit. Thus it would be an interest- 
ing experimental challenge to look for signature of the 
AGFP at intermediate temperatures. 

Before we conclude, it is worth mentioning that the 
geometry studied in our paper is of direct interest for the 
production of non-local entangled electron pairs propa- 
gating in two different wires. These electron pairs are 
produced by Cooper pair breaking via crossed Andreev 
processes when the superconductor is biased with re- 
spect to the wires comprising the junction. One can 
ask if electron-electron interaction in the wires actually 
leads to enhancement of entangled electron pair pro- 
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duction via the crossed Andreev processes. For exam- 
ple, we have observed that for the NSN junction with 
r — 0,t — 0,rA 0,tA ^ 0, interaction can lead to 
enhancement of the crossed Andreev amplitude (i^) un- 
der RG flow. This implies that for the case when the 
superconductor is biased with respect to the wires, inter- 
electron interactions enhance the production of non-local 
entangled pairs, for which the amplitude tA is relevant, 
as compared to local entangled pairs for which the am- 
plitude rA would be relevant. This is consistent with 
the results of Recher and Loss^ who also argued that 
it is energetically more favourable for the two entangled 
electrons of the Cooper pair to go into different wires, 
rather than the same wire. Finally the RG flow leads to 
a fixed point with = 1 where the system becomes a 
perfect entangler. A more general case would be when 
r ^ O.t ^ O.TA 0,tA ^ 0- To study this case, one 
can start from the two wires SFPS-matrix and study the 
RG flow of tA for an S-matrix which is in the close vicin- 
ity of this fixed point. The result of this study is shown 
in Fig. [HI We show that starting from the short-distance 



cut-off d the RG flow initially leads to enhancement of 
Ia which will lead to an enhancement in the production 
of non-local entangled pairs. Hence, these studies estab- 
lish the fact electron-electron interactions the wires can 
actually lead to an enhancement of non-local entangled 
electron pair production. 
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APPENDIX A 

Here we give the RG equations for nine independent 
parameters in case of a FSN junction. 
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APPENDIX B 



Here we give the RG equations for six independent pa- 
rameters in case of a symmetric 3 wire NSN junction. 
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dt 
ll 
dt' 
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dl 

dt^A 
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|[r*(ri + t^+t'^ + tl + t'l) - r{l - \r\^)] - a'[r|r^p + r\{ttA + t't'^)] 
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